First VLTI/MIDI observations of a Be star: a Ara 
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Abstract. We present the first VLTI/MIDI observations of the Be star a Ara, showing a nearly unresolved 
circumstellar disk in the N band. The interferometric measurements made use of the UTl and UT3 telescopes. The 
projected baselines were 102 and 74 meters with position angles of 7° and 55° , respectively. These measurements 
put an upper limit to the envelope size in the N band under the Uniform disk approximation of </>max ~ 4± 1.5 mas, 
corresponding to 14 i?*, assuming i?*=4.8R0 and the Hipparcos distance of 74 pc. 

On the other hand the disk density must be large enough to produce the observed strong Balmer line emission. 
In order to est imate the possible circumstellar and stellar parameters we have used the SIMECA code developed 
by Stee l)l995h and Stee & Bittar l)200lf) . Optical spectra taken with the echelle instrument Heros and the 
ESO-50cm telescope, as well as infrared ones from the 1.6m Brazilian telescope have been used together with the 
MIDI spectra and visibilities. These observations put complementary constraints on the density and geometry 
of a Ara circumstellar disk. We discuss on the potential truncation of the disk by a companion and we present 
spectroscopic indications of a periodic perturbation of some Balmer lines. 

Key words. Techniques: high angular resolution - Techniques: interferometric - Stars: emission- line. Be - Stars: 
winds, outflows - Stars: individual [a Ara) - Stars: circumstellar matter 



1. Introduction 

Be stars are hot stars exhibiting the so-called Be phe- 
nomenon, i.e. Balmer lines in emission and infrared ex- 
cess. Optical/infrared and ultraviolet observations of Be 
stars have been widely interpreted as evidence for two 
quite distinct regions in the circumstellar environment of 
these objects: a rotating, dense equatorial region, called 
"disk" in the following, and a diluted polar region which 
expands with veloc ities that may reach 2000 kms~^ (e.g. 
Marlborough Il986|) . 

The thermal infrared domain is an important spec- 
tral region where the transition between optically thin 
and thick disk occurs. The optically thick disk emission 



is growing from 7 to 15 /im. It reaches about half of the 
total continuum flux at 8 /im and dominates the Spectral 
Energy Distribution (SED) almost completely at 15-20 /Ltm 
(Gehrz et al. .1974) . The flux from the central star at 
these wavelengths is 10 to 50 times fainter compared to 
A = 2 /im whereas the ratio Fcnv/F* reaches a factor 4- 
10. In particular for nearly edge-on disks, the flux from 
the optically thick disk is very sensitive to the inclination 
angle, since it is proportional to the emitting surface. 
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* Based on observations made with the Very Large Telescope 
Interferometer at Paranal Observatory 



Observations at 20 /im by Gehrz et al. demon- 
strated that the mid-IR excess is due to free-free emis- 
sion fro m the hydrogen envelope. The ad-hoc model from 
Waters lll986l) has been successful to explain near- and 
far-IR observations and is coherent with polari zation 
data (Cote & Waters ^9 87. Wate rs fc M arlborough |l992L 
Dougherty et al. ll994l Cote et al. ll99(j) . In this model the 
IR excess originates from a disk with an opening angle 



tl and a density distribution p(r) oc r~ seen at an ob- 
server's angle i. From IRAS data of four Be stars Waters 
found that the far-IR slope of the SED is not strongly 
influenced by the opening angle, as long a.s 9 + i < 90°. 
Moreover, he found that n—2A gives a good agreement for 
X Oph, S Cen and (p Per if the disk is truncated to 6.5 R*. 

Interferometry in the visible has also been used to 
study the circumstellar e nviron ment of Be stars ( Thorn 
et al.[l9.86, Mourard et al. ll989l Quirrenbach et al. ll994l: 
Stee et al.|l995|j). Quirrenbach et al. l)l997|) gave an upper 
limit on the opening angle of about 20° and using spec- 
tral Differential Interferometry (DI) structures within Be 
disks were monitored with a hi gh sp atial resoluti on dur - 
ing several years (Vakih et al. Il998l Berio et al. \l99^ . 
Nevertheless, the disk extension at IR wavelengths is still 
subject to an active debate since various authors give quite 
different exten sions as a function of wavelength. For in- 
stance Waters l|l98fil) found a typical IR extension of 16 i?^ - 
from 12, 25 and 60 /im observations, Gehrz et al. (1974) 
found 8i?Vf at 2.3 and 19.5 /xm, Dougherty et al. (jl994j) 
a disk size lar ger th an 20 i?* in the near-IR, and finally 
Stee & Bittar Ij200ll) obtained that the Br7 emission line 
and the nearby continuum originates from a very extended 
region with a typical size of 40 i?*. 

In order to solve this issue, we have used the 
VLTI/MIDI interferometer to obtain the first IR measure- 
ments in the N spectral band of a Be star. The selected 
target, aAra (HD 158 427, HR651G, B3Ve) is one of the 
closest Be stars with an estimated distance of 74 ± 6 pc 
derived from the Hipparcos parallax, and is well known as 
emission line star s i nce t he discovery of its II/3 emission 
by Pickering l)l896l ll'896|) . aAra is a 2.9 magnitude star 
in the Johnson V-Band and its IRAS flux at 12 micron is 
12.7 Jy. Its color excesses E{V - L) and E{V - 12 /xm) 
are respectively 1.8 and 2.23, being among the highest of 
its class. The projected rotational velocit y v sin i has been 
estim ated t o be 250-300 km s^^ (Yudin l200]l Chauville 
et al. l200l|) . An intrinsic linear polarization of Pl«0.6% 
with a Position A ngle (P A) of 172° h as also been detected 
(McLean & Clarke ll979l Yudin l200ll) . Since the disk orien- 
tation is expected to be perpendicular to this direction, we 
expe ct the disk ma ior-axis to be aroun d PA » 82° (Wood 
et al. ll996alll996al Quirrenbach et al. ll997(l . 

The spectral type definition of a Ara is constant in 
the literature oscillating between B2Ve and B3Ve. We 
ado pted B3 Ve following the latest study of Chauville et 
al. I|2001[) . The typical stellar radius and effective temper- 
ature for this spectral type are 4.8 Rq and Tcff = 18 000 
respectively. 

By using IRAS observations of aAra, Waters l|l986t) 
estimated the outer disk radius is estimated to be about 
7 R*, i.e. 2.1 mas at a distance of 74 pc. Stee 12003 ) pre- 
dicts the visibility of a Ara at 2 /im to be lower than 0.2 
with a 60m baseline, i.e. fully resolved. Using the same 
model parameters, a significant visibility loss should still 
be detectable in the N band. This paper will show that at 
10 /im a Ara is in fact unresolved which, combined with 
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spectroscopic data, puts some constraints on density, in- 
clination angle, disk orientation and distance of a Ara. 

The paper is organized as follows. In Section |2] we 
present the interferometric MIDI observations and the 
H-band spectroscopy. Section 01 briefly describes the 
SIMECA code and the envelope parameters used for the 
modeling of a Ara. The modeling strategy is exposed in 
details. In Section^jwe present our theoretical results and 
discuss the possibility for a close but unseen companion to 
truncate the disk making a Ara unresolved in the N band. 
SectionElpresents some spectroscopic variations of the II/3 
line profile that may be another evidence for this hypo- 
thetical companion. Finally, Section|Hldraws the main con- 
clusions from these first IR interferometric measurements 
of a Be star. 

2. Observations 

2.1. Interferometric data 

We use the VLTI / MI DI I nterferometer described by 
Leinert et al. l|2003al l2003b|) which combine the mid-IR 
light (N band) from the two VLT Unit Telescopes Antu 
(UTl) and Melipal (UT3). The observations of a Ara were 
done during the nights of June, 16th and 17th 2003 and 
are presented in Tabled 



Un the loth ot June, a Ara was observed with a 1U2 m 
baseline and a PA=7° whereas June, 17 the projected 
baseline was 79 m and PA=55°. 

The observing sequ ence, described extensively in 
Przygodda et al. l)2003^ is summarized hereafter. The 
chopping mode (f=2Hz, angle —90°) is used to point and 
visualize the star. The number of frames recorded by im- 
age are generally 2000, the exposure time is by default 
4 ms in order to avoid the quick background saturation. 
If the result of the acquisition sequence is not satisfac- 
tory, this procedure is re-executed. Then, the MIDI beam 
combiner, the wide slit (0'.'6x2"), and the NaCl prism are 
inserted to disperse the light and search for the fringes by 
moving the VLTI delay lines. The resulting spectra have 
a resolution A/ A A ~30. 

When searching for the fringe signal, the delay line of 
the VLTI is moved, while the MIDI internal piezo-driven 
delay line is performing additional scans of the fringe pat- 
tern. Once the fringes are found another file is recorded 
while MIDI is tracking the fringes by its own, i.e. by per- 
forming a real time estimation of the OPD based on the 
data continuously recorded. The correction is send to the 
VLTI delay lines at a rate of about IHz. Finally, the pho- 
tometric data are recorded in two more files. 

In the first file, only one shutter is opened, correspond- 
ing to the calibration of the flux from the first telescope, 
here UTl, and the flux is then divided by the MIDI beam 
splitter and falls on two different regions of the detector. 
The total flux is determined separately by chopping be- 
tween the object and an empty region of the sky, then 
the source flux is computed by subtraction. Once it has 
been done, the same procedure is carried out again for the 
beam arriving from UTS. 

We performed the spectral reduction and fringe cal- 
ibration using a software developed at the Max-Planck 
Institut fiir Astronomic in Heidelberg, writte n in IDL^. 
The reader is referred to Leinert et al. l)2004j) for an ex- 
tensive discussion on MIDI data reduction and error han- 
dling. 

The photometric datasets used for the calibration of 
the contrast of the dispersed fringes are also used for the 
spectrum flux calibration. The frames obtained on the tar- 
get and the ones obtained on the sky are averaged. Then, 
the averaged sky frame is subtracted from the averaged 

^ Package available at this address: 

|http://www.mpia-hd.mpg.de/MIDISOFT/| 

Table 2. Calibrator stars parameters, from the MIDI cal- 
ibrators catalogue, (Van Boekel et al. 2005). 
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Fig. 1. a Ara flux measured by the VLTI/MIDI (8-13 /im) 
instrument with ISO data (2.5-11.5 /im) plotted as solid 
line compared with the theoretical SED of the star alone 
(dashed line) and the star-|-envelope emission (dotted line) 
assuming a distance of 105 pc. Note the IR excess starting 
around 5 /im, i.e. log A=0.7. 

target frame. In the resulting integration, the spectral axis 
is oriented along the horizontal detector axis. For each de- 
tector column, a Gaussian function is fitted to the peak. 
The position of the spectrum in all illuminated columns, 
as function of the column number, is fitted by a quadratic 
polynomial. In a similar way, the width of the spectrum is 
fitted by a linear function. This procedure is carried out 
on both photometric datasets, corresponding to telescopes 
UTl and UTS, respectively. Both fits are then averaged 
and used to create a mask consisting of a Gaussian func- 
tion with the average position and width of the spectra for 
each column. Finally, a weighting function is computed to 
create a mask for the extraction. 

The photometric masks are also used to extract the 
fringes information. Each frame of the fringe data is re- 
duced to a one-dimensional spectrum by multiplying it 
with the mask; integrating in the direction perpendicu- 
lar to the spectral dispersion; and finally subtracting the 
two output channels of the beam combiner from each 
other. The spectra from each scan with the piezo- mounted 
mirrors are collected and Fourier-transformed from opti- 
cal path difference (OPD) to frequency space. The fringe 
amplitude at each wavelength is then obtained from the 
power spectrum. We typically co-added the signal by bin 
of four pixels in the dispersion direction to improve the 
signal-to-noise ratio. 

The scans where fringes were actually detected are 
selected based on the white-light fringe amplitude, i.e. 
the fringe amplitude that is seen after integrating over 
all usable wavelengths. The histogram of all white-light 
fringe amplitudes within a fringe track dataset usually 
shows a small peak near zero, and a broad peak at higher 
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Fig. 2. Ha and H/3 line profiles observed in mid-1999 (solid lines) compared with the best SIMECA fits (dotted lines). 



amplitudes. We interactively set a threshold just below 
this broad peak, and average the fringe amplitudes of all 
scans with a white-light fringe amplitude higher than this 
threshold to give the final fringe amplitude as a function of 
wavelength. The raw visibility is obtained by dividing the 
fringe amplitude by the total flux. The calibrated visibility 
is obtained by dividing the raw visibility of an object by 
that of a calibrator. The photometrically calibrated flux 
creating the fringes is called correlated flux. 

The interferometric and photometric data are pre- 
sented in Table 1 and the parameters of the interferometric 
calibrators are given in Table 2. 

The flux calibration has been performed with the data 
of June, 16 only, a Aql was considered as a good calibra- 
tor for the absolute flux extraction and we have used the 
ISO observations of this star as a flux reference (but has 
not been used as interferometric calibrator). An indepen- 
dent calibration is performed for the individual spectra 
from each part of the detector and each telescopes. The 
airmass is also corrected independently for each of them 
by using the observations of one calibrator for two differ- 
ent times. Fig. n shows the VLTI/MIDI Spectral Energy 
Distribution (SED) with the ISO data and the theoretical 
SED obtained with the SIMECA code (see hereafter). 

The errors on the observed visibilities are mostly sys- 
tematic. The statistical signal-to-noise ratio on the white 
light fringe amplitudes is 5-10 at minimum and much bet- 
ter after adding up the several hundred scans taken per 
interferometric measurement. Most of the time, the dom- 
inant error source is the fluctuations of the photometry 
which occur between the fringe recording and the photom- 
etry recording from the individual telescopes which are 
separated in time by 2 to 10 minutes. This affect the fringe 
signal for object and calibrator measurements. Comparing 
the raw visibilities observed for different calibrator stars 
during one night, the standard deviation of these values 
under good conditions amounts to ± 5-10 % (relative) 



at the red and blue end of the spectrum while for adverse 
conditions these numbers have to be multiplied by a factor 
of 2-2.5. These fluctuations are much larger for instance 
than the systematic error expected from the diameters of 
the c alibrators and are mostly achromatic (Leinert et al. 
l2004l) . Yet, the instrumental visibility curve of MIDI is 
very stable in shape, i.e. the slope of the curve is typically 
varying by less than 3-5%. The visibilities observed on cal- 
ibrators ("instrumental visibility") rises from about 0.4 at 
8 [ivn to about 0.7 at 13 /xm, rather repeatable from night 
to night (i.e. it means that, in good atmospheric condi- 
tions, the absolute level is dependant on the photometric 
variations and hence the stability of the atmospheric con- 
ditions but the shape of the visibility curves is mostly an 
instrumental parameter only weakly affected by the atmo- 
sphere) . 

The data reduction described previously does not con- 
sider the variations of the photometry during the 1- 
2min needed to record the data and the visibility esti- 
mation is usually based on the mean photometric flux. 
Moreover, under good atmospheric conditions, the statis- 
tics of the correlated flux per individual scan follow 
roughly a Gaussian curve and the mean of this correlated 
flux is used for the visibility estimation. The main origin of 
the variability of the correlated fiux is the high frequency 
photometric fluctuations, the atmospheric piston and the 
quality of the fringe tracking (i.e. the instrumental piston). 

During the night of June, 17 2003 thin cirrus were 
passing through the telescope beams. The seeing was quite 
stable, below 0.5" during all the night but the standard de- 
viation of the flux from the target pointed by the (visible) 
seeing monitor^ was oscillating between 0.02 and 0.05 be- 
tween 3h and 8h (UT time) , indicating the passage of thin 
cirrus. These fluctuations seem correlated with an acceler- 

^ This information has been extracted from the ESO 
Ambient conditions database of Paranal observatory: 
http: //ar chive. eso.org/i 



atlon 01 the atmosphere turbulence with a coherence time 
in visible decreasing from to=:5.5 ms at 3h UT to tq—2 ms 
at 7h30 UT. By comparison, the photometric fluctuations 
of the seeing monitor calibrator during the 16th of June 
are below 0.01 during the entire night although the turbu- 
lence was also quite rapid with tq—3 ms. The variability 
of the atmosphere during the June, 17 night was such that 
a careful inspection of the stability of the photometry and 
the correlated flux recorded has been necessary. The pas- 
sages of the thin cirrus were easily recognized by a rapid 
drop of the flux lasting between 0.5 and 2s during the pho- 
tometric record and by a simultaneous photometric and 
correlated fluxes drop in the tracking files. The frames 
and the scans concerned were excluded from the data re- 
duction process: this represent about 60% of the photom- 
etry and 50% of the scans, i.e. about 800 frames were left 
from each photometric files and about 100 scans were used 
for the visibility estimation. Of course, this treatment, ap- 
plied to a Ara and its calibrators is a first order correction 
and the standard deviation of the instrumental visibility 
curves of the calibrators observed between 2h UT and 9h 
UT is still large compared to good observing night, on the 
order of 18% whereas the standard deviation of the previ- 
ous night is only 8%. Fortunately, the longest baseline of 
102m has been used during the best night. The error bars 
reported in Fig. Eland Fig.Elreflect the standard deviation 
computed by this mean. 



2.2. Spectroscopic data 

For an investigation of the variability of a Ara this study 
makes use of a set of 42 echelle spectra, taken during 69 
nights in May to July 1999. The data were secured with 
the Heros instrument, attached to the ESO-50cm tele- 
scope. The resolving power is A A/ A — 20 000, recording 
the regions from 350 to 565 nm and from 585 to 860 nm in 
a blue and red channel simultaneously. The spectra have 
a typical signal-to-noise ratio of 170 (Fig|21). 

Next to these data, complementary spectra have also 
been recorded quasi-simultaneously with the VLTI run in 
order to know whether a Ara had shown IR emission lines. 
We have observed a Ara in the J2 band (1. 2283-1. 2937/im) 
using the 1.6 m Perkin- Elmer telescope and Coude spec- 
trograph(with R=10 000) at the Observatorio do Pico dos 
Dias, Laboratorio Nacional de Astrofisica (LNA), Itajuba, 
Brasil. In total, five stellar spectra were recorded with 
the Camara Infravermelho (CamlV) detector, August 13, 
2003. Each spectrum corresponds to a different position of 
the star along the spatial axis of the entrance slit, allow- 
ing to use the median of the five two-dimensional frames, 
from which the average bias frame was subtracted, as a sky 
background. After dividing the stellar and sky background 
frames with the average flat field, the sky was subtracted 
from two-dimensional stellar spectra. The extraction of 
one-dimensional spectra, as well as wavelength calibration, 
using an Ar-Ne calibration lamp, have been performed 
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Spectral type 

Mass 
Radius 
Luminosity 



B3 Ve 
18 000 K 
9.6 Mq 
4.8 Rq 
2.2 10^L(7 



0_ 



with standard IRAF^ packages. Finally, the continuum 
was normalized to unity. 

The observed Pa/3 line proflle is shown in Figure |31 
The emission intensity was I / Ic ~ 2.3, while the violet- 
to-red peak height ratio V/R was 1.17. Even if the spec- 
troscopic data were not recorded fully simultaneously with 
the interferometric ones, available observations of typical 
timescales of the spectral changes of a Ara, such as the 
above described Heros data, justify the assumption that 
the circumstellar environment of a Ara has not changed 
signiflcantly between June and August 2003. This Pa/3 line 
proflle, combined with the interferometric data, constrains 
the physical parameters of the circumstellar environment 
of a Ara, as it will be shown in the following sections. 

3. The SIMECA code and the envelope 

parameters used for the modeling of a Ara 

In order to constrain the physical parameters of the 
circumstellar environment of a Ara, we have used the 
SIMECA code. Th is cod e, describe d in p revious paper 
(see Stee & Araujo Il994t Stee et al. Il995l: Stee & Bittar 
l200lh . has been developped to model the environment of 
active hot stars. SIMECA computes line profiles, Spectral 
Energy Distributions (SEDs) and intensity maps, which 
can directly be compared to high angular resolution ob- 
servations. The envelope is supposed to be axi-symmetric 
with respect to the rotational axis. No meridional circu- 
lation is allowed. We also assume that the physics of the 
polar regions is well repr esente d by a CAK type stellar 
wind model (Castor et al. Il975|) and the solutions for all 
stellar latitudes are obtained by introducing a parameter- 
ized model which is constrained by the spectroscopic and 
interferometric data. The inner equatorial region is dom- 
inated by rotation, therefore being quasi Keplerian. The 
ionization-excitation equations are solved for an envelope 
modeled in a 410*90*71 cube. In order to reproduce the 
Pa/3 profile we have computed the line corresponding to 
the transitions between the 5-3 atomic levels. The popu- 
lations of the atomic levels are strongly altered from their 
Local Thermodynamic Equilibrium (LTE) values by non- 
LTE conditions in the wind. For computation, we start 
with the LTE populations for each level, and then com- 
pute the escape probability of each transition in the wind. 



IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of 
Universities for Research in Astronomy (AURA), Inc., under 
cooperative agreement with the National Science Foundation. 



obtaining up-dated populations. By using these popula- 
tions as input values for the next step, we iterate until 
convergence. The basic equatio ns of t he SIMECA code 
are given in detail by Stee et al. (1995). 

To take into account the photosphcric absorption line, 
we assume the underlying star to be a normal B3 Ve star 
with Toff = 18 000 K and R = 4.8Ro and synthesize the 
photosphcric line profil es using the SYN SPEC code by 
Hubeny (Hubeny Il988l: Hubeny & Lanz The re- 

sulting line profile is broadened by solid-body rotation and 
might be further altered by absorption in the part of the 
envelope in the line of sight towards the stellar disk. 

The stellar parameters (Table I^J are important mainly 
to define the distance and the luminosity of the source 
(Fig-CJ E^nd represent the first step of the iterative process 
leading to a model of the envelope as described below. In 
the scope of this model, we have computed the Ha, H/3 
(FigEl and Pa/3 (Fig|2|) line profiles. 

The critical constraint of this study is based on the fact 
that the Ha and H/3 lines (FigI2l, as well as the the Pa/3 
line (Fig. are strongly in emission, whereas from Fig. 
the visibilities are clearly around unity, i.e. the envelope 
in the N band is nearly unresolved. Hence, the circum- 
stellar density must be large enough to produce Balmer 
lines in emission while at the same time this density must 
be low enough and/or the circumstellar environment must 
be well delimited for keeping the envelope unresolved. The 
emissive lines of many Be stars like a Ara exhibit a char- 
acteristic double peaked shape and the 'violet' and 'red' 
peaks are denoted V and R. This shape provide stringent 
constrains on the wind parameters and particularly on the 
inclination of the system. 

In order to obtain matching parameters derived from 
both line profiles and interferometric measurements we 
have adopted the following strategy: 

1. The density po at the base of the wind is fixed, having 
the greatest influence on the emission line intensity. 

2. A range of inclination angles is chosen in order to 
match the V and R peak separation, V and R peak 
intensities and line wings with the observed profile. In 
this case, the range is 40-50°. 

3. In the next step, the ml and m2 parameters are varied, 
where ml describes the variation of the mass flux from 
the pole to the equator according to: 



m - (/-pole + [(0cq. - ^/-polc) Sin"l(0)] 



(1) 



where is the stellar colatitude. The parameter m2 
describes the variation of the terminal velocity v^o as 
a function of the stellar latitude: 

voo{0) = Ko(pole) + [voo(eq.) - Voo(pole)] sin"^^^^) (2) 

where Voo {pole) and Voa (eg) are respectively the polar 
and equatorial terminal velocities. 
4. The ratio between the equatorial and polar mass flux: 



Cl= 



pole 



(3) 



IS typ ically between lU and lU (Lamers and Waters 
Il987 ) , having a large influence on the wings of the line 
profile. 

5. Following the above procedure, the Ha and H/3 line 
profiles are fitted simultaneously to obtain the general 
envelope parameters. 

6. With the parameters reproducing well the 1999 ap- 
pearance of Ha and H/3, the modelled Pa/3 line profile 
has a larger intensity than actually observed in 2003 
(Fig. 13). In the following, we assumed that the only 
parameter that might have changed between 1999 and 
2003 is the density at the base of the envelope. In fact, 
with a base density decreased by about 25 %, a bet- 
ter agreement between the observed and modeled Pa/3 
profiles is obtained (Fig. OJ. 

7. The comparison between the theoretical SED we ob- 
tained and the MIDI and ISO data (Fig. ^ is used 
to constrain the distance of a Ara. The adopted a dis- 
tance is based on the hypothesis on the star parameters 
and the fitted parameters of the circumstellar environ- 
ment of a Ara which settle the infrared excess and the 
visual reddening (see discussion Sec. 15.21) . 

8. In a final step, we compute the intensity maps in spec- 
tral channels of 0.2 /im. Using the distance estimated 
in the previous step, we obtain theoretical visibilities 
which are compared to the VLTI/MIDI data. We found 
that the critical parameters to fit the observed visibil- 
ities are the distance to the star, the inclination angle, 
the ml value and the density at the base of the wind 
Po- 

During the above described procedure, more than 500 
models were computed using the SIVAM II Alpha ES45 
EV68 workstation, build on 3 servers with each four 1 
GHz CPU and 24 GB of RAM at the Observatoire de 
la Cote d'Azur. To compute one model it takes about 
20 min of CPU time. 

The variability of the circumstellar environment of 
aAra is traced by the variations of the Ha line equiva- 
lent width (EW) reported in the literature by Dachs et 
al. il990); Hanuschik et al. l)l995|) and this work. From 
Fig. 21 it is evident that we cannot straightforwardly use 
parameters derived from the 1999 Ha and H/3 line profiles 
to model the Pa/3 line profile observed in 2003, i.e. close to 
our interferometric run. It is most likely that some physi- 
cal parameters have changed between these two epochs. 



4. Results and Discussion 

4.1. Theoretical line profiles and SED 

The Ha line observed 1999 shows strong emission of 
I/Ic w 3.7, which is well reproduced by SIMECA with 
the stellar parameters given in Table |31 and the envelope 
parameters from Table 0] (Fig. 121. 

With a wind base density decreased by 25%, produc- 
ing a good agreement for the Pa/3 line profile (Fig. OJ, 




Fig. 3. Left: Observed Pa/3 line profile (niid-2003, solid line) compared to the modeled one based on the parameters 
determined by the fit of the 1999 Ha and H/3 line profiles (dotted line). Right: Observed Pa/3 line profile (solid line) 
and the best fit from the SIMECA code (dotted line). To compensate for the effects of non-simultaneous observations, 
the density of the disk was allowed to vary in a small range. This fit has been obtained with a disk density lower by 
25% compared to the best model based on 1999 Ha and H/3 lines. 




Fig. 4. Variation of the Ha equivalent width (E W) wit h 
time. The data were taken from Dachs et al. l|l99Cll) : 
crosses; Hanuschik et al. (,1995) : stars and this work: dia- 
mond. 

the model gives an Ha line profile with a I/Ic ratio «3 
(Fig. ISJ. Unfortunately, we do not have Ha line profiles 
taken contemporary to the VLTI observations that could 
confirm such decrease. 

Even though Fig. 13 shows that the Pa/3 line profile 
modelled with SIMECA is in good agreement with the ob- 
served one, the modelled V/R is less than unity, whereas 
the observed one is V/R w 1.17. SIMECA is based on 
a radiative wind model for which the gas is outflowing 



Fig. 5. Ha line profile obtained in April 1999 with the 
HERDS instrument at la Silla (Chile) and the SIMECA 
fit (dotted line) for the base density used to fit the 2003 
Pa/3 profile. The strength of the line is comparable of the 
one reported in 1982 or 1987, for instance see Fig. 2 by 
Hanuschik et al. lll995.1 



from the star with a polar and equatorial velocity of re- 
spectively 2000 and 170 kms"'^. In order to obtain a 
V/R > 1, it would be necessary to introduce an asym- 
metry like a global one-armed oscillation, as seen by Berio 
et al. l|l999l) and Vakili et al ( 199a) for the 7 Cas and 
C Tau equatorial disks. Nevertheless, the global agreement 
of observation and model both in profile shape and inten- 



sity indicates that, under the general assumptions ot the 
SIMECA model, the global kinematics and density dis- 
tribution within the envelope are well reproduced by the 
assumed model parameters. 

The computed SED between 2 and 10 /Ltm is shown in 
Fig.n together with the data recorded by VLTI/MIDI (8- 
13.5 fim) and the ISO satellite (2.5-11 /zm). The observed 
SED and the modelled curve match for a source distance 
of 105 pc. This distance is significantly larger than the 
one obtained from the Hipparcos satellite (see discussion 
in Sect.EJ. In Fig.^the continuum emitted by the central 
star is approximated as a Teff=18 000 K black body radi- 
ation plotted as a dashed line. The total, i.e. free-free and 
free-bound emission from envelope -I- star, is indicated by 
the dotted line. Clearly, the infrared excess produced by 
the circumstellar envelope is mandatory in order to fit the 
observed data (solid line). 

4.2. Theoretical visibilities 

The parameters obtained by the above modeling of the 
emission line profiles, and the density required to fit the 
2003 Pa/3 line, can now be used to compute the expected 
visibility curves for a distance of 105 pc, as estimated from 
the fit of the SED. The resulting visibility curves are plot- 
ted in Fig. IHl clearly showing that the modeled envelope 
should be well resolved at 79 m and 102 m (mean visi- 
bility Vmean ~ 0.63) whereas the VLTI/MIDI data with- 
out any doubt have hardly resolved the target {Vmean ~ 
1). However, the observed visibility curves in Fig. may 
indicate that the envelope is more resolved for the longer 
wavelengths. This effect, if true, is more or less reproduced 
by our model and is more obvious for the 102 m baseline 
where the theoretical visibility is 0.62 at 8 fim and 0.57 at 

13 IITlH. 

Nevertheless, we did not succeed to obtain model pa- 
rameters reproducing both the Balmer lines and the visi- 
bility curves at the same time: At a distance of 105 pc the 
model gives an envelope which is clearly resolved, if the 
constraint that the Pa/3 line is in emission with I/Ic ~ 2.3 
is met. One possibility would be to increase the distance 
of the star up to 122 pc, a distance suggested by Cohen 
et al. l|l998l) . However, we consider this value as too high 
by comparison to the Hipparcos distance (see also the dis- 
cussion in Sec|3J). 

The modeled envelope is more or less spherical since 
the input parameter ml = 0.3 corresponds to an opening 
angle of 160° (see Stee 1200^ . producing only a negligi- 
ble visibility difference for the envelope seen parallel or 
perpendicular to the polar axis. Nevertheless, this small 
flattening may produce the small intrinsic linear polariza- 
tion of Pl«0.6% measur ed by McLean & Clarke (197^ 
and Yudin et al. l)l998|) . who found a polarization an- 
gle PA=172°. The inner disk orientation is, therefore, ex- 
pected to be perpendicular to this direction at 82°. This 
means that the longer baseline was unfortunately oriented 
at about 70° from the expected major axis of the disk. 



Ihe diagram ot the a Ara circumstellar envelope projected 
onto the sky plane (with an arbitrary oblateness) and the 
VLTI baseline positions is given in Fig. El 

Considering these difficulties, we have computed a 
model with a truncated disk. The radius where the trun- 
cation occurs was set to 22 R^,, which was derived using a 
distance of 105 pc and the fit of the Pa/3 line. (Fig.[7I. The 
theoretical visibilities in the Ha, H/3, Pa/3 and Br7 lines 
are plotted in Fig. IHlfor the two scenarios. For both mod- 
els it is clear that at shorter wavelengths, with a 102 m 
baseline, the envelope of a Ara would be well resolved in 
the Pa/3 and Br7 lines {V « 0.2-0.3) even for a truncated 
disk {V « 0.35-0.45). In Ha and H/3, V is 0.6 and 0.85, 
respectively, for the full disk and V is 0.65 and 0.94 for the 
truncated one. At these wavelengths, the emission comes 
from the inner part of the disk, which remains unresolved 
at a distance of 105 pc. Consequently, such a truncation 
of the disk would have only a little impact on the Ha and 
H/3 emitting regions, whereas the IR emission, originating 
from more extended regions, is more affected. 

The fact that a Ara is nearly unresolved at 8 /um with 
a 102 m baseline gives an upper limit for the diameter 
of the emitting envelope of (j>rnax= 4±1.5 mas, i.e. 14 Rj, 
(for d —74 pc) or 20 R* (for d — 105 pc) assuming an 
uniform disk distribution for the star -f disk brightness dis- 
tribution. If a gaussian distribution is assumed, the best 
fit is obtained with a FWHM=2.5mas although the qual- 
ity if the fit is poorer than for the uniform disk hypothe- 
sis. Such a radial extension is in good agreement with the 
one obtained by Waters et al. l|l987h from IRAS Far-IR 
data, with (j)disk '--^ 7 R^. Indeed, the hypothesis that the 
star-|-envelope flux distribution is unchanged from 8/im 
to 13.5/im is a crude approximation and these estimated 
extensions must be taken as indicative only. It must be 
stressed out that the slopes of the visibility curves pro- 
vided by the model are in good agreement with the data 
showing that the flux is more concentrated at 8/im than 
at 13^m. Moreover, the angular diameter estimations are 
mostly based on the visibilities taken with the longest pro- 
jected baseline at PA=7° to the East, i.e. more or less 
along the predicted polar direction (see Fig.EJ. 

In Table 01 we present the results based on the best fit 
of the a Ara August 2003 Pa/3 line profile and visibilities 
in the N band. 



5. Discussion 

The characteristics of the disk around a Ara seem to have 
been fundamentally stable over the past century. Only the 
overall strength of the emission, as traced e.g. by the Ha 
EW exhibit large variations (see Fig. ^J. However, the 
compilation of spectroscopic observations and the consec- 
utive study of this disk by MIDI have revealed a number 
of interesting and partly unsuspected characteristics. 
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Fig. 6. VLTI/MIDI data for the 79 m (left) and 102 m (right) baselines and theoretical visibilities ft-om SIMECA 
(dashed line) for the distance of 105 pc estimated from our SED fit (see Fig.^). The error bars are equivalent to one 
sigma. 
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Fig. 7. VLTI/MIDI data for the 79 m (left) and 102 m (right) basehnes and theoretical visibilities fr-om SIMECA 
(dashed line) for a truncated disk at 22 R^. by a possible companion. The error bars are equivalent to one sigma. 



5.1. Variability and suspected binarity 

The Hydrogen recombination line profiles of a Ara ex- 
hibits (quasi?)-periodic short-term V/R variations with 
a timescale of 45 days (Mennickent & Vogt, 1991). The 
timescale for these variations measured in the HEROS 
dataset, however, is better accounted for with a larger pe- 
riod of 70-90 days (see below). Superimposed to these vari- 
ations are changes of the emission line equivalent widths. 
For instance, the strength of the Ha emission varies from 
about 3.2 to more than 4.5 in units of the local continuum 
in only a few weeks. The period of variations could be the 
same as the V/R period, but this cannot be shown from 
our data and it also has not yet been noticed by other 



investigators. The authors studying the V/R variations 
concentrated most of their efforts on the shape of the line 
profile rather than on the absolute strength. Nevertheless, 
some interesting cases of rapid V/R changes, correlated 
with strong variations of Balmer line equivalent widths 
are reported by Mennickent (1989; 1991). 

Long-term V/R variations are common in Be stars, be- 
ing attributed to a one-armed oscillation in the Keplerian 
disk. Rapid V/R variations, however, have not been ex- 
tensively discussed in the literature. Mennickent (1991) 
associates fast quasi-periodic variations (P<100d) to the 
smallest disks, i.e. the disks with the lowest equivalent 
widths in Ha. Okazaki (1997) found that for early-type 
Be stars around BO radiative effects can explain the con- 




Fig. 8. Theoretical visibilities for the Ha, H/3, Pa/3 and Br7 lines for baseline orientations corresponding to the ones 
observed on June 16, i.e. PA= 7° onto the sky plane, for both scenarios discussed. On the left side we show model 
data for a disk with the parameters which fits the 2003 Pa/3 line profile, but not in agreement with the VLTI/MIDI 
visibilities (Fig. 0). On the right side model data for a disk truncated at 22 R^, e.g. by a putative companion, is 
presented, being in agreement with the observed visibilities (Fig.Ej). 



parameter /result 


value 


uncertainty 


Distance 


105 pc 


10 pc 


Inclination angle i 


45 


5 


Photospheric density (pphot) 


1.2 10-i^g cm"^ 


linked to CI 


Photospheric expansion velocity 


0.07 km s"^ 


0.01 


Equatorial rotation velocity 


300 km s"^ 


20 


Equatorial terminal velocity 


170 km s"^ 


20 


Polar terminal velocity 


2000 km s"^ 


500 


Polar mass flux 


1.7 IO-^Mq year"^ sr"^ 


0.5 10"'' 


ml 


0.3 


0.1 


w2 


0.45 


0.05 


CI 


30 


linked to Pphot 


Mass of the disk 


2.3 10"^°Mq 




Mass loss 


6.0 10"^Mq year"^ 





Table 4. Best model parameters for the a Ara circumstellar environment obtained by fitting the visibilities and the 
2003 Pa/3 fine profile with SIMECA 



finement of this one-armed oscillations with typical peri- 
ods of about 10 yr. Nevertheless, this model would have 
difficulties to explain more rapid V/R variations. Another 
possibility for such variable asymmetry is the presence 
of a companion, orbiting at distances larger than the en- 
velope extension. Then, these quasi-periodic line profile 
variations should be somehow related to the orbital pa- 
rameters of the putative secondary. 

Nevertheless, with respect to the photospheric absorp- 
tion lines a Ara is unusually stable among earlier type 
Be stars. In a sample of 27 early Be stars observed with 
Herds observations, only two stars, including a Ara, did 
not show line any detectable profile variability (Rivinius 
et al. 2003). 

However, this is quite different for the circumstellar 
emission lines. Next to slow, possibly cyclic long-term 



changes typical for Be stars, there are changes on a shorter 
timescale. In particular, the peak-height-ratio of the violet 
and red peaks of the higher Balmer lines have made one 
full cycle of low amplitude (about 0.95 to 1.05) during the 
69 days of the Herds observations (see Fig. El- At the 
same time, the radial velocity of the emission component 
of the Balmer lines was changing in a cyclic way as well. 
For instance. Fig. 1111 shows this effect for the central ab- 
sorption of the H/3-emission, but the same is seen if the 
radial velocity of the emission peaks is traced, or if other 
Balmer lines are investigated. Unfortunately, the available 
data are not good enough to claim such a radial velocity 
change in the weaker emission lines or in the photospheric 
absorption lines. The uncertainty of measuring the radial 
velocity of such lines is much higher than the amplitude 
seen in Fig. Mennickent et al. (1991) have observed 
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Fig. 9. Diagram of a Ara circumstellar envelope pro- 
jected onto the sky plane (with arbitrary oblateness) based 
on polarization (P) me asure ments (PA=172°) o btaine d 
by McLean & Clarke l|l979l) and Yudin et al. (Il998^. 
The actually observed VLTI baseline positions are shown 
for June, 16: PA = 7°, B==102m and June, 17: PA ^ 
55°, B=79m following the polarization (P) measurements 
(PA=172°) obtain ed by McLean & Clarke ((1979,1 and 
Yudin et al. (Il998l) . 

similar changes and gave a cycle length of 0.13 years, or 
45 days. However, their Fig. 2 would also support a cycle 
length of 70 to 80 days, as present in the Heros data of 
1999. 

Such behavior is known from a few other Be stars as 
well. In the cases investigated thoroughly, it was shown to 
be linked to binarity, although the mechanisms are not yet 
known (see for instance Koubsky et al. 1997). A search for 
any spectral contribution of such a hypothetical compan- 
ion in the phase-binned Heros spectra of a Ara did not 
return a positive result. 

Nevertheless, since a companion as a potential origin 
for an outer truncation of the disk is important in light of 
the above results, the properties of such a system are esti- 
mated from the following: Using the stellar parameters of 
Tabled i.e. a mass of 10 M©, a 70 day period would give 
a radius of about 154 R0, assuming a circular orbit of a 
companion with negligible mass. With i?* = 4.8 R©, this 
corresponds to about 32 stellar radii, which is in agree- 
ment with the estimate based on the MIDI/ VLTI data for 
disk truncated at 25 R*, i.e. somewhat smaller than the 
companion orbit. 

Such an orbital dimension can also be estimated by 
integrating the radial velocity amplitude detected in II/3 
of the order of 12-16 km s~^. Assuming a circular orbit, a 
period of 69 days and testing two inclination angles oii = 
30° and 60°, the mass of the companion would range from 
1.4 Mq (F2-4V) and 2.8 M© (A2-4V). The corresponding 
Roche lobe radii RRoche are roughly between 15 and 20 R^, 



althougli wc point out that the radiai velocity curve shown 
in Fig^] would suggest an eccentricity of about 0.2-0.3, 
affecting the previous estimation. 

However, the radial velocities derived from emission 
lines cannot be taken at face value. Depending on the 
systemic properties, the radial velocities measured in the 
central inversion of the circumstellar emission follow the 
actual radial velocity curves of the stars only loosely. The 
above orbit size and period, again for the circular case 
and negligible mass of the companion, result in an or- 
bital velocity of about 18kms^^, which is in the same or- 
der of magnitude as the amplitude measured in the emis- 
sion lines. In some well- investigated binaries like 59 Cyg 
(Harmanec et al. 2002, Rivinius & Stefi 2000, P-29d) or 
(/) Per (Hummel & Stefl, 2001, Stefl et al. 2000, P-127d), 
the emission of the Balmer lines is indeed in phase with 
the orbital motion. 

Until now, no evidence has been reported for a pos- 
sible companion around a Ara, but it should be kept in 
mind that it is difficult to search for low mass compan- 
ions around bright stars, such as Be stars (see for instance 
the case of 7 Cas, Harmanec et al. 2000; Miroshnichenko, 
Bjorkman, & Krugov 2002, period P~204d). The X-ray 
flux from a Ara is not peculiar and similar to that of nor- 
mal B stars ( Ly=3 .75 lO^^erg s~^, assuming d=122pc'*, 
Cohen et al. Il998l) . The proposed truncation may also 
help to explain the stringent non detection of radio emis- 
sion from a Ara at a level of 0.1 mJy at 3.5cm and 6.3cm 
(Steele et al. 1998). 

Using the NPOI interferometer, Tycner et al. l|2004^ 
have recently studied the disk geometry of the Be star 
C Tau, which is also a well-investigated spectroscopic bi- 
nary (P'^133d, K~10kms~^). They measured the disk 
extension quite accurately to be well within the Roche ra- 
dius. This suggests also that this disk may be truncated. 

5.2. Stellar parameters and distance estimation 

The distance used to model the star deviate from the 
Hipparcos distance by several sigma. The error bar of the 
Hipparcos measurement is 6 pc, and even taking into ac- 
count a putative correlated error of th e ord er of 1 mas^ 
(Narayanan & Gould, 1999, Pan et al. \2004 . the largest 
distance within 1 sigma could be 87 pc. The distance of 
105 pc comes straightforwardly from chosen stellar pa- 
rameters, i.e. the absolute magnitude as estimated from 
effective temperature and radius. We notice for instance 
that Fabregat & Reglero (1990) obtain a distance of 98 pc 
from direct ubvy photometry. Zorec & Briot (1991) obtain 
85 pc, and note that their estimates of Be star distances 

^ This implies a luminosity Lx=1.02 lO^^erg s~^, assuming 
the Hipparcos distance d=74pc, which is still in the range of 
lumin osities encountered for this spectral type by Cohen et al. 
il998h 

^ This correlated error is realistic since Narayanan & Gould, 
(1999) report spatially correlated errors of Hipparcos up to 
2 mas in the direction of the Pleiades and Hyades clusters. 
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Fig. 10. Variation of the height ratio of violet and red 
emission peaks of the H/3 Una during the Heros observa- 
tions 1999 
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Fig. 11. Radial velocity changes of the central depression 
of H/3 during the Heros observations 1999 
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Fig. 12. Average of the lUE-spectra of aAra (dotted) 
compared to a theoretical spectrum with Tcff = 18 000 K 
and log (7 = 4.0 (solid line). 




seem systematically lower by 20 pc compared to normal B 
stars observed in clusters, suggesting that their distance 
is a lower limit. Also comparing the flux-calibrated lUE 
archival spectra to Kurucz ATLAS9-data gives a distance 
of about 105 pc (see Fig. IT^ . assuming R^, — 4.8 R0 and 
neglectable extinction in the UV. The latter is probably 
justified since there is no interstellar extinction towards 
such a close star (Dachs et al., 1988), and the circumstel- 
lar effects for a non-shell star (see discussion below) should 
not affect the UV-regime. We note also that the param- 
eters of a Ara and its hydrogen environment provided by 
the SED fits of Dachs et al. (1988) are close to the ones 
presented in this paper, in particular the estimated cir- 
cumstellar reddening. 

The available spectral classifications of the star are all 
around B3Ve, so that it can be assumed that the effec- 
tive temperature is well constrained at about 18 000 K. 
Instead, the main source of uncertainty in the model may 
be the radius, which arbitrarily has been chosen to be 



Fig. 13. Schematic view of the a Ara circumstellar envi- 
ronment as used in this work (see Table EJ. 

4.8 Rq, based on the statistical value given by Tokunaga 
in Allen's "Astrophysical Quantities" (2000). 

However, aAra has a high wsini of about 270 km s^^ 
(Dachs et al. 1990, Yudin 2001, Chauville et al. 2001). 
With such a large projected rotational velocity, the hy- 
pothesis of a uniform disk radius is questionable and the 
radius used is most probably underestimated as demon- 
strated by the recent observations of aEri (Achernar, 
Domiciano de Souza et al. 2003). This implies an increase 
of the illuminating area and strengthens the difficulty to 
put the star at the Hipparcos distance. Without consid- 
ering any reddening, and keeping the Hipparcos distance, 
the radius of aAra would be unrealistically low (below 
3.5 i?©) or the photosphere unrealistically cold {Tcs of the 
order of 15 000 K). In Tab. [S] we show selected parameters 
of the Be stars a Ara, a Eri and the 'normal' (i.e. non Be) 
star TyUMa which are all classified as B3V. Comparing 



these values, it becomes evident that, ii their Hipparcos 
distances are correct, the striking visual flux differences 
between these three stars cannot be directly related to a 
distance difference only. 

The reddening of Be stars does not only contain an 
interstellar term, which is zero for a Ara, but also has a 
circumstellar contribution. This is a model dependant pa- 
rameter, explaining the well known difficulty to calibrate 
the Be stars distances independently from Hipparcos mea- 
surements. 

For instance, Cohen et al. l)l998(l have apparently to- 
tally neglected any circumstellar reddening^ and derived 
a distance of 122 pc. This difficulty also affects the present 
model and an underestimation of the circumstellar extinc- 
tion gives a distance too large. On the other hand, the 
measured B — Y for a Ara does not indicate a circumstel- 
lar environment dense enough in the line of sight to explain 
the difference, since this should be linked to emission line 
properties which are not observed. 

A large amount of extinction and reddening of the cen- 
tral star would have to be due to an edge-on disk, and 
such an inclination would lead to a large and deep self- 
absorption in the Balmer emission lines, the so-called shell 
appearance, a Ara, however, has never been seen in a shell 
phase. Instead, the inclination angle is relatively well con- 
strained by the shape of the Balmer line to a value of 
about 45°. 

An inclination of about 45° implies an equatorial ro- 
tational velocity of the order of 380kms~^, since wsinz = 
270kms~^. Using the parameters in Tabled this value 
is at about 75% of the critical velocity. We point out 
that the equatorial velocity in the model was considered 
a free parameter and the best model provides an equato- 
rial value of 300 km s^^, i.e. 60% of the critical velocity. 
With such a high rotational velocity, the von Zeipel effect 

® We note also that their estimated distances for other Be 
stars like & Cen seem also systematically overestimated com- 
pared to Hipparcos distances. In contrary their estimated dis- 
tance of oEri is only 27pc based on a wrong assumption of 
the radius of this star (Re5=Rpoi=4.1R0) which is indeed very 
difficult to constrain without any interferometric observations. 

Table 5. Stellar parameters of three close B3V stars, 
a Ara, a Eri and r\ UMa. 





a Ara 


a Eri 


jyUMa 


V 


2.9 


0.5 


1.7 


Dist.' 


74.3±6pc 


44.1±1.4pc 


30.9±0.7pc 


V sin i 


270 kms"^ 


225 km s"^ 


150 kms~^ 


B-V 


-0.17 


-0.17 


-0.17 


V-K" 


0.4 


-0.4 


-0.5 


lv,Gq 


4.8 R0 


12 ± 0.4 Ro 


4.8 R© 


■r> Hi 


4.8 Ro 


7.7 ± 0.2 R© 


4.8 Ro 


Hipparcos distance. 



" From 2MASS pliotomctry . 

4.8 is the statistical value for tliis a B3V star, tire radius for a Eri is 
from Domiciano de Souza et al. 2003. 



IS becoming important, and thus the continuum emission 
must be strongly latitude dependant (see e.g. Townsend et 
al., 2004, his Fig. 3). This may also confuse the line pro- 
file diagnostic and affect the angle determination. Such 
an effect is not included in the present model. A value 
of the inclination closer to 90°, compatible with the large 
t;sini observed, may be related to a stronger circumstel- 
lar reddening from the equatorial environment and a lower 
integrated stellar ffux since the equatorial regions of the 
star are cooler and fainter than the polar ones. 

Definitely, the Hipparcos distance, if reliable, is a very 
constraining parameter for any model of the circumstellar 
environment of a Ara. More interferometric observations 
are needed, firstly to constrain the angular size and shape 
of the hydrogen emission lines, and secondly to evaluate 
the angular size and shape of the underlying star. The 
forthcoming VLTI/ AMBER instrument will be able to ad- 
dress the first point by studying the fringe properties in 
Pa/? and Br7. However, baselines on the order of 150-200 m 
will be mandatory to resolve the central star, becoming 
available only with the 1.8m Auxiliary Telescopes. Their 
first light is foreseen for 2005. 

6. Conclusion 

The first IR interferometric measurements of a Ara sug- 
gest that the size of the circumstellar envir onment is 
smaller than what was predicted by Stee l)2003|) . The fact 
that a Ara remains unresolved, but at the same time 
exhibits a strong Balmer emission line puts very strong 
constraints on the parameters of its circumstellar disk. 
Independently of the physical model we used, these mea- 
surements put an upper limit of the envelope size in the 
N band of (?!)max= 4 mas, i.e. 14 R^, if the star is at 74 pc 
according to the Hipparcos parallax or 20 if the star is 
at 105 pc as suggested by our model, assuming a spher- 
ical uniform disk for the star-t-envelope brightness distri- 
bution. The suspected presence of a low mass companion 
could help to understand the limit of the disk extension 
as determined by MIDI. A schematic view of the a Ara 
circumstellar disk, as it was modeled in this work with the 
SIMECA-code, is presented Fig. ESI 
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